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C^ ' Abstract. We study the survival of a single diffusing lamb on the positive half line 

C^ . in the presence of N diffusing lions that all start at the same position L to the right of 

' ^ I the lamb and a haven at x = 0. If the lamb reaches this haven before meeting any lion, 

Y^ • the lamb survives. We investigate the survival probability of the lamb, Sn{x^ L), as a 

• i-H , function of N and the respective initial positions of the lamb and the lions, x and L. 

^^' We determine Sn{x,L) analytically for the special cases of iV = 1 and N — >■ c». For 

<~i ■ large but finite N, we determine the unusual asymptotic form whose leading behavior is 

, Mh | Sn{z) ^ N~^ , with z = xjL. Simulations of the capture process very slowly converge 

to this asymptotic prediction as N reaches 10^°". 
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1. Introduction 

We investigate the one-dimensional diffusive capture process in which a marlced particle 
— a "lamb" — diffuses on the positive half line a; > in the presence of N independently 
diffusing predators — "lions" — that are all initially at L > x. If the lamb meets any 
lion, the lamb is killed. Additionally, the origin is a haven for the lamb. If the lamb 
reaches the haven before meeting any of the lions, then the lamb survives. We are 
interested in the survival probability of the lamb as a function of the starting positions 
of the two species, as well as on the number of lions. 

This model is a natural counterpoint to the well-studied capture process of a 
single diffusing lamb in the presence of N independent, diffusing lions on the infinite 
line [H [21 |3l m |5]. In the most interesting situation where the lions are all on one 
side of the lamb, the survival probability SNit) of the lamb asymptotically decays as 
a power-law in time, Si\f{t) ~ t^'^^, with the exponent P^ exhibiting a non-trivial 
dependence on the number of lions A^ and also on the diffusivities of each animal. For 
simplicity, the case where the diffusivities of all animals are the same (and set to one) 
is normally considered. The initial positions of the lamb and the lions are irrelevant in 
this asymptotic behavior. 

For this diffusive capture on the infinite line, the exponent /3jv is known exactly only 
for iV = 1 and AT = 2: /^^ = 1 and /^s = f [D El S H El El [7] . The latter result shows 
that even though the two lions are independent, their effect on the capture process is 
not, since t~^^ > {t^^^) ■ For the case A^ = 3, a mapping to an equivalent electrostatic 
problem leads to the accurate estimate /^g = 0.91342 ± 0.00008 |8j. For A^ > 3, the 
value of Pn has been estimated with moderate accuracy only for a few values of A^ [1] ; 
however, it is known that (3^ > 1 for A^ > 3, so that the average lifetime of the lamb is 
finite [9]. Because P^ grows more slowly than linearly with A^, each additional lion has 
a progressively weaker influence on the capture process. As A^ — > oo, both asymptotic 
and rigorous arguments give /^at — )■ j\nN [21 [3l H] . Parenthetically, the capture process 
with lions sited on both sides of the lamb is much more efficient than in the one-sided 
system. For A^ lions with approximately equal numbers of them on either side of the 
lamb, the lamb survival probability asymptotically decays as t~^^ , with 7^^ growing 
linearly with A^ for large A^. 
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Figure 1. Space-time trajectories of a lamb and two lions with a haven at a; = 
when: (a) the lamb survives and (b) the lamb is eaten. 

In this work, we incorporate the new feature of a haven at a; = and ask whether 
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the lamb can reach the haven before meeting any of the hons. If the haven is reached, 
we say that the lamb survives (Fig. [1]). Our goal is to determine how the ultimate 
survival probability Sn{x,L) depends on the initial positions of the lamb and all the 
lions, X and L > x, respectively, as well as on the number of lions. As we shall see, 
the survival probabihty depends on z = j^ rather than on x and L separately and thus 
we write the ultimate survival probability as Sn{z). Our main result is that Sn{z) has 
an unusual form whose leading behavior is Sn{z) ~ N~^ , but this behavior does not 
become apparent until A^ becomes of the order of 10^°*^. 

We begin by solving the simplest and exactly- soluble case of one lion in Sect. [2j 
We also outline the formal solution to the problem for any number of lions. In Sect. [3] 
we treat the extreme case where the number of lions is infinite, so that the lion that 
is closest to the lamb moves ballistically. We then investigate arbitrary A^ in Sect. HI 
When A^ is large, we can replace the A^ lions by a single "closest lion" that moves 
deterministically. We develop approximation schemes to estimate S^iz) in this large- 
A^ limit. We also present numerical results for the survival probability in Sect. O A 
straightforward simulation of the random-walk motion of the particles is prohibitively 
slow when A^ is large, and we present two alternative approaches that are considerably 
more efficient and allow us to probe the survival probability in the regime where A^ is 
extremely large — of the order of 10^°°. Finally, in Sect. [6l we summarize and also 
discuss some natural and intriguing extensions of the model. 

2. Exact Analysis 

2.1. One Lion 

As a preliminary, we can readily solve the case of one lamb at x = Xi and one lion 
aX L = X2 > Xi for the general situation where the diffusivities of the two species are 
distinct — Di for the lamb and D2 for the lion. We compute the survival probability 
that the lamb reaches the haven at x = before being eaten by the lion, S{xi,X2), 
by mapping the coordinates of the lamb and the lion on the line to diffusion in a two- 
dimensional wedge, from which the survival probability follows easily. 

It is convenient to transform from the coordinates (xi,X2) to yi = X\l ^fTTx and 
?/2 = X2/ \/Ty2- In the yi-y2 plane, the motions of the lamb and lion on the half line 
can be viewed as the isotropic diffusion of a fictitious composite particle with unit 
diffusivity [TJ [10]. If yi reaches zero while the condition yi < ?/2 is always satisfied, the 
lamb survives (Fig. [2]). Conversely, if y\\fD\ = y2y/D2 at some time (corresponding to 
Xi = X2) while yi always remains positive, then the lamb has been eaten by the lion 
before the haven is reached. 

In the yi-y2 plane, the initial position of the composite particle is 

Xi X2 
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Figure 2. Mapping the diffusion of a lamb and a lion on the half line a; > to isotropic 
diffusion in a wedge of opening angle 6. 



corresponding to the polar angle 
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a = tan 



(1) 



The allowed region for the composite particle is a wedge of opening angle 

9 = tan-i VD2/D1 . (2) 

We want the probability S{xi,X2) that the composite particle first hits the line yi = 
(corresponding to the lamb reaching the haven) without hitting the line yi^/lDi = 
y2y/D^- This probability satisfies the Laplace equation [10] 










dXi- UX2 

for X2 > xi, with boundary conditions S{xi = 0,X2) = 1 and S{xi,X2 = xi) = 0. Clearly 
the solution is a function that linearly interpolates between and 1 in the angular 
direction, so that the ultimate survival probability is [TOl [11] 

a - (7r/2 - 6) 
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Figure 3. Ultimate survival probability S{xi^X2) as a function of X\lx2 for various 
values of the lamb and lion diffusivities, D\ and £'2, respectively. 

As is obvious from Fig. |3l the closer that the lamb starts to the haven the more 
likely it is to survive. Moreover, as can be inferred from Fig. [21 the best strategy for the 
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lamb to survive for a given initial condition is to diffuse quickly. As the diffusivity of the 
lamb Di increases, the wedge angle in Fig. |2] approaches ^ while the starting position 
of the fictitious particle in the plane moves close to the 1/2 axis, i.e., closer to the haven. 
Finally notice that in the limit D2 ^ (stationary lion), the survival probability decays 
linearly with X\lx2- 

As a byproduct of the wedge mapping, we can immediately determine the 
probability that the lamb is still diffusing — that is, the lamb has not yet reached 
the haven and has not yet been eaten by the lion. This situation corresponds to the 
fictitious particle having not yet reached either of the sides of an infinite wedge defined 
by xi = and x\ = X2. In the isotropic yi-y2 coordinates, this wedge has opening angle 
6 (Fig. [2]), and the survival probability asymptotically decays as t~'^l'^^ . In particular, 
when D\ = D2, then 6 = j (see Eq. (J2])), and the survival probability asymptotically 
decays as t~^. 

2.2. Formal Solution for General N 

The reasoning given above can be readily generalized to map the problem of a diffusing 
lamb in the presence of N diffusing lions to a single diffusing fictitious particle in A^ + 1 
dimensions, with boundary conditions that reflect the lamb reaching the haven or being 
eaten by a lion. For simplicity, we set the diffusivities of the lamb and the lions to 
one. We first discuss the case of two lions; the generalization to any number of lions is 
immediate. 

Suppose that the lamb is initially at Xi > and that the two lions are initially at 
X2 = X3 > Xi. The lamb survives if it reaches x = without meeting either of the lions 
on the way to x = 0. We now map the diffusion of the three interacting particles on 
the positive half line to the isotropic diffusion of a composite particle at {xi,X2,X3) in 
three dimensions, with constraints that correspond to the interactions in the lamb-lion 
system. By this mapping, the allowed region for the composite particle is defined by 
xi > 0, corresponding to the lamb not yet reaching the refuge, as well as by xi < X2 and 
Xi < X3, corresponding to the lamb not yet eaten by either of the lions. This defines a 
wedge-shaped region that are delineated by three planar sides that is known as a Weyl 
chamber [12] . 

The survival of the lamb corresponds to the composite particle first hitting the 
plane xi = of the Weyl chamber without hitting either of the planes xi = X2 and 
Xi = X3. By the equivalence between first-passage and electrostatics fTO], this survival 
probability of the lamb coincides with the electrostatic potential $(xi,X2,X3) at the 
initial point of the composite particle, with the boundary conditions $ = 1 on the 
plane Xi = 0, and $ = on the planes Xi = X2 and Xi = X3. This same mapping 
works for any number of lions and constitutes the formal solution. Unfortunately, the 
analytical solution to this potential problem does not seem tractable for more than one 
lion (i.e., three or more particles), although some extreme value electrostatic properties 
have recently been exactly solved for the three-particle problem |1] 
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3. Infinite Number of Lions 

When the number of hons is infinite, the lion closest to the haven — the closest lion — 
would reach the haven at an infinitesimal time. However, it is instructive to consider the 
related problem in which each lion undergoes a nearest-neighbor random walk. In this 
case, the position of the last lion inexorably moves one lattice spacing to the left in each 
time step. For this system, we determine the ultimate survival probability Soo{x,L) by 
writing the backward Kolmogorov equation ^U\ for 5*00(3;, L) and then applying scaling 
to solve this equation. The result should correspond to that obtained for diffusing lions 
in the limit of very large N. 

To write the backward equation, we consider the evolution of the system over a 
small time interval [0, At] during which the lamb moves to x + ri{0)At and the boundary 
moves to L — vAt, where v is the boundary velocity. That is, the position of the lamb 
x{t) evolves by the Langevin equation dx/dt = i]{t), where ri{t) is Gaussian white noise 
with zero mean, {ri{t)) = 0, and correlation {'ri{t)ri{t')) = 2 D 6{t — t'). We now view 
the new positions of the lamb and the boundary after the time interval At as the initial 
conditions for the subsequent evolution. Thus S{x,L) = {S{x + rj{fS)At,L — vAt)), 
where the average is over the initial noise ^^(O). Expanding the right-hand side of this 
recursion to lowest non-vanishing order in each variable and using the properties of 
delta-correlated noise, we obtain the backward equation 
d^S dS 

for < a; < L, with the boundary conditions 5(0, L) = 1 and S{L, L) = 0. To solve 
this equation we make the scaling ansatz S{x, L) = f{y) (with y = x/\/L) to give the 
ordinary differential equation for /: 

r + ^yf' = o, (5) 

subject to the boundary conditions /(O) = 1 and f{vL) = 0; here the prime denotes 
differentiation with respect to y. 

Integrating and applying the boundary conditions gives 

/(,) = 1 - -^f (^)/^^) . (6) 

^^ erf(vWlD) ^ ^ 

In the limit L — )■ 00, this expression reduces to 

f{y) -^ erfc(y v/^V4^) = eTic{z^/vL/AD) , (7) 

with z = x/L. The primary feature of this result is that the lamb survival probability 
is non-zero only within a thin boundary layer where the starting position satisfies 
X <ti \/ADL/v. Outside this layer the lamb is almost surely eaten by one of the lions. 

4. Asymptotics for Large A^ 

The capture process also simplifies when the number of lions A^ is finite but large, 
because the position of the closest lion becomes progressively more deterministic as A^ 
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increases, even though each individual hon undergoes independent Brownian motion. 
Thus we only need to consider the ultimate survival of the lamb in the presence of a 
single effective predator p^ — the closest lion — that moves systematically towards 
the lamb (Fig. H]). We now exploit this physical picture to give a heuristic argument for 
the ultimate survival probability of the lamb. 




Figure 4. Schematic space-time representation of the diffusion of the lamb at x and 
the position L{t) of the closest lion when the number of lions A^ ^ 1. 



When all the lions start at L > x, the average number of lions at x is 



n{x,t) 



-{x-LY/ADt 



V4jrDi 

We estimate the location of the closest lion, L{t), by demanding that n{L{t),t) = 1. 
This criterion gives |1] 

L{t) = L-\/At, (8) 



where 



A = AD\yiM\ 1 



llnlnM 



2 1„M+-j ('^ 

and M = N/y/^. Thus to lowest order, A ^ ADlnN. At a critical time t* = L^/A 
the closest lion has reached the haven at x = and the capture process is necessarily 
finished — either the lamb has been killed or it has reached the haven. Notice that 
although N must be large for the closest lion to move deterministically, A^ cannot be 
too large. As discussed in the previous section, if each lion undergoes a nearest-neighbor 
random walk, the closest lion moves deterministically to the left with speed v = 1 when 
A^ is sufficiently large. For Eq. ([8]) to be valid, we therefore require (to lowest order) that 
y/ADt InN < vt, or N < exp^v'^t/AD). Using v = 1 and Z) = | for a nearest-neighbor 
random walk, the last lion moves deterministically as v4JDtlnN only when t > 2 In A^. 
For t < 21nA^, the last lion moves with constant unit speed toward the lamb. 

We now crudely estimate the ultimate survival probability of the lamb as the total 
probability flux that reaches x = up to time t* in the semi-infinite system without 
any additional constraints. This integrated flux represents an upper bound for the 
survival probability for large A^ because this estimate includes lamb trajectories that 
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could intersect the trajectory of the last lion and then reach the haven. For a diffusing 
particle that starts at x, the flux to an absorbing boundary at the origin at time t is: [10] 

Consequently, the probability Sn for the lamb to get trapped at the origin up to time 
t* (corresponding to the lamb reaching the haven and surviving) satisfies the bound 



X 



Sn< e-'/'^^^^t =eTic{zVh[N). (10) 

Jo V47rDt3 



Here we have used the substitution u = x/y/ADt to transform to a Gaussian integral, 

x_ 



as well as the lowest-order approximation t* = L'^/{ADlnN) and z — ^ 



From the asymptotic form erfc(x) ~ e ^ /(vvnr^), we thus obtain an upper bound 
for the ultimate survival probability that has the unusual functional form 

Sn < ^ N-^' [ln{N'')]-'/\ (11) 



Consistent with basic intuition, Sn is a decreasing function of A^ and also decreases 
as z — > 1 with A^ fixed. It should be emphasized that Eq. ( TTTl) applies in the limit of 
zV^nN ^ 1, which is extremely hard to achieve by direct simulation. For example, if 



4 



the lamb starts halfway between the haven and the lions (z = |), then for A^ = 10 
the argument of the complementary error function is zy/\nN f« 1.52; for A^ = 10^^, 
zy/hiN ^ 3.03. Conversely to reach zy/hiN = 10 requires A^ = e^°° ~ 10^'''^. Notice also 
that Eq. ( TTOj) matches the survival probability given by Eq. (^^ for a ballistically-moving 
boundary when A^ reaches a critical value for which the completion time t* = L'^/A also 
equals L/v. 

More rigorously, we should also incorporate the absorbing boundary condition at 
L{t), corresponding to the lamb getting eaten by the closest lion. This problem of a 
fixed absorbing boundary at x = and a moving absorbing boundary at x = L — VM 
does not seem readily soluble, however. Instead, we investigate a related model in 
which the boundary motion mimics that of the closest lion, but is engineered to be 
soluble. As we shall show, the ultimate survival probability for this alternative problem 
has a qualitatively similar dependence on system parameters as Eq. ( fTTl) . Consider 



the toy model in which the closest lion coordinate is Ltoy(t) = y/L'^ — Bt (compared to 
L(t) = L — y/At, with A = ADlnN, for A^ ^ 1 diffusing lions). These two boundaries 
satisfy the inequality Ltoy(t) > L{t) and both reach the origin at the same time when 
B = A. Thus the toy model remains an upper bound for the true survival probability. 
It is again convenient to treat the evolution of the system in the two-dimensional 
space {x,L). Let S{x,L) be the probability that the lamb successfully reaches the 
haven, where x and L denote the initial positions of the lamb and the boundary 
respectively. Following the same approach as in Sec. [3l we write the backward equation 
for S{x,L). In a small time interval [0, At] the lamb moves to x -|- ri{0)At, where rj{t) 
is Gaussian white noise with zero mean, and the boundary moves to L — {B/2L)At, 
where B/2L is the boundary speed. The survival probabihty now satisfies S{x, L) = 
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{S{x + ri{0)At, L — {B/2L)At)), and expanding the right-hand side to lowest order gives 
the backward equation 

dx^ 2LdL ^ ^ ' 

for < X < L, with the boundary conditions S{0,L) = 1 and S{L,L) = 0. To solve 
( !T2|l we make the scaling ansatz S{x, L) = f{y), with y = 72;, where 7 = ^JB/{2D) and 
z = ^, and find that the scaling function satisfies 

r+z//' = o, (13) 

for < y < 1, with the boundary conditions /(O) = 1 and /(7) = 0; here the 
prime denotes differentiation with respect to y. Integrating once gives / oc e~^ '^, 
and integrating again gives 

f{z) = 1 - erf (7z/v^)/erf (7/v^) , (14) 

where the constants are determined by the boundary conditions. Substituting in 
7 = -^ B /{2D) and B = AD In A^ gives the asymptotic behavior 

f{z) ^ [N-'" - iV-i] ~ N-'' . (15) 

This upper bound has the same asymptotic behavior as ( TTTl) and suggests that the 
heuristic approach should be quite accurate. 

5. Simulations 

We now present simulation results for the lamb-lion-haven system. While a direct 
simulation is simple to code, it becomes prohibitively slow when N is large. We have 
therefore developed two complimentary approaches to determine the survival probability 
in the large- A^ limit. 

5.1. Probability Propagation 

Probability propagation is well-suited for probing the case of iV ^ 1, where we replace 
the position of the closest lion by a deterministic absorbing boundary, L(t), that moves 
according to Eq. (|8]). Here, the constant A can be chosen as the mean or most probable 
position of the closest lion or any other reasonable positional metric. We choose to 
set A = ADlnN, which is the leading behavior in Eq. (|8]). The omission of higher- 
order corrections, which slightly decrease A, lead to a more slowly-moving boundary 
and a correspondingly slightly larger survival probability. Thus probability propagation 
should provide a lower bound to the true survival probability. 

Let P{x, t) be the probability that the lamb is at x at time t. At each time step, 
the probability in the interior region 2 < x < \_L{t)\ — 1 propagates according to 
P{x,t -|- 1) = ^P{x — l,t) -|- |-P(x -|- l,t); here [i^(t)J is the largest integer less than 
L{t). At the edge sites P(l,t+1) = |P(2,t) and P([L(t)J, t + 1) = iP([L(t)J -l,t). 
Probability elements that reach either a; = or [i^(t)J -|- 1 do not propagate further and 
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remain in place. Probability propagation continues until L{t) reaches x = 0. The total 
probability at a; = at this termination time gives the survival probability of the lamb. 
We used probability propagation to obtain Sn{z) for N up to 10^°°. We used 
quadruple precision variables to ensure accuracy of the probability values throughout 
the propagation. The initial value of L was chosen to be the smallest such that finite-size 
effects were imperceptible — this ranged from L = 1000 for small N to L = 30, 000 for 
the largest A^ values. 

5.2. Event-Driven Simulation 

A naive simulation is simply to move every lion and the lamb by ±1 at each time step, an 
approach which is prohibitively slow for large A^. However, there is no need to simulate 
every single random-walk step, particularly if the lamb is far from both the haven and 
nearest lion. This motivates using an event-driven simulation, in which we propagate 
all particles over a time that corresponds to a finite fraction of the time needed for a 
reaction to actually occur — either the lamb reaching the haven or getting eaten by the 
closest lion. 

Let y be the minimum of the distances between the lamb and the nearest lion, and 
between the lamb and the haven. We could move every particle according to a binomial 
distribution oiy/2 — 1 steps because there is no possibility that the lamb meets any of the 
lions or reaches the haven during this update. However, this approach is unnecessarily 
stringent because each particle moves a typical distance that is only of the order of y/y. 
Thus we increment the number of steps by m, where 

(y'/2Y iiy>Y, 
m = { 16 

\y/2 iiy<Y, 

and move every particle according to a binomial distribution of m steps. Note that 
these update rules match at the crossover separation y = Y . After each such update, 
we check if the lamb has reached or crossed over the position of the haven or that of 
any lion, in which case the simulation is finished. 

For y < y , 771 = I and the lamb cannot reach either the haven or any lion during the 
update; this part of the simulation is exact. For y >Y , there is a non-zero probability 
that the lamb trajectory could cross the haven or a lion trajectory and then cross back 
during the update. However, by choosing Y appropriately, the probability of error due 
to such crossing trajectories can be made vanishingly small. We found that Y = lb 
gave an excellent compromise between accuracy and efficiency. We also checked that 
simulations results with the update rule f lT6|) are essentially identical with exact results 
that arise by choosing y = oo in the update rule ( 1T6l) . 

5.3. Results 

In Fig. |5]^a) we show the dependence of the ultimate survival probability versus scaled 
initial position z = j^ for A^ up to 256, 000, with 10^ realizations for each data point, from 
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event-driven simulations. For A^ ^ 1, these survival probabilities gradually converge, 
as A^ increases, to a limiting curve that corresponds to the system where the last lion 
moves ballistically. The lions are all initially at L = 100 and we verified that the survival 
probability depends only on the ratio x/L, without any explicit finite-L dependence. 
This independence on L emerges when L > 100 and thus we focus on the smallest 
system (L = 100) where finite-size effects are negligible. 
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Figure 5. (a) Dependence of the survival probability on z for representative values 
of N from event driven simulations. Dashed curve is analytic solution for ballistic lion 
motion from Eq. ([6|) where v = 1, D = 1/2, and L = 100. (b) Survival probability 
versus number of lions N for three representative z values. Curves give the analytic 
prediction S = N^^ and the symbols represent probability propagation results. 



We also examined the dependence of Sn on A^ for fixed z to test the asymptotic 
power-law behavior S ~ A^~^ of Eq. (ITT|) . Our analytical prediction matches the 
simulation quite well for z < 0.5 (Fig. [5])(b). However, a small but slowly growing 
discrepancy arises as z is increased beyond 0.5. The source of this discrepancy is that 
the heuristic derivation of Sec. H] ignores the existence the absorbing boundary caused 
by the last lion. When z approaches 1, the lamb starts sufficiently close to the last lion 
that the assumption of ignoring the boundary caused by the last lion is no longer valid. 

Finally, we compare our two simulation approaches with each other and with our 
our heuristic prediction Sn ^ erfc(z-\/ln A^) from Eq. ( ITOl) . By construction, the event- 
driven simulation is more accurate because it explicitly follows the stochastic motion 
of the lamb and the lions. Our heuristic prediction flTOj) provides an upper bound for 
large A^, but this regime is not feasible to simulate with the event-driven algorithm. 
Conversely, the probability propagation simulation can be implemented for arbitrarily 
large A^ but suffers from systematic error because it assumes the closest lion position to 
be deterministic. 

Figure [6t^a) illustrates the convergence of the simulation results to Eq. flTOj) where 
the difference between the simulated value of Sn{z) and erfc(2;\/ln A^) is plotted as 
a function of z for representative A^ values. We quantify this difference by A = 
{Ai^ — As)/Ag_, where A^, = j Sjq{z)dz, with Sn{z) = erfc(zA/ln A^), is the area beneath 
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Figure 6. (a) Difference between simulated survival probability and Sn{z) = 
erfc{zy/\nN) as a function of z for various N. Open symbols correspond to probability 
propagation while filled symbols correspond to event-driven simulations. Circles 
correspond to iV = 8, 000, triangles to TV = 32, 000, and squares to iV = 128, 000. (b) 
Relative area difference A versus 1/lniV for probability propagation (o) and event- 
driven (♦) simulations. 



the analytic survival curve and similarly for the area beneath the simulated curve. 
Figure EJ^b) shows that A — )■ as iV — )■ oo for the probability propagation algorithm. 
A similar, but not identical convergences arises in the event-driven simulation, but the 
method cannot reach the large-iV regime. These results provide strong evidence that 
the survival probability is indeed given by 5 — ?■ erfc(2;\/lniV) as N -^ oo. 

6. Outlook 



The presence of a haven adds an intriguing element to the classic capture process of 
a single lamb in the presence of A^ diffusing lions. Now the basic question is whether 
the lamb can reach safety at the haven before it is eaten by one of the lions. We 
investigated the dependence of the ultimate survival probability of the lamb, Sn{x, L), 
on the number of lions A^ and also on the initial positions of the lamb (x) and the 
lions (all at L, for simplicity). By a rough heuristic argument, we found that Sn has 
the asymptotic behavior Sjy ^ erfc(2A/ln A^), and this function has the unusual leading 
behavior Sn ~ N~^ , where z = x/L. It is remarkable that a simplistic approach gives 
such an unusual and rich result. However, the approach to this asymptotic regime is 
extremely slow and it is necessary to simulate a system that corresponds to A^ of the 
order of lO^"*^ lions before the asymptotic behavior becomes apparent. 

It is natural to ask about the properties of the ultimate survival probability in 
higher dimensions. For diffusive capture in an unbounded system, the case of one 
dimension is the most interesting. However, the presence of a haven now makes the 
higher- dimensional problem nontrivial. For example, in two dimensions, a natural 
setting would be a diffusing prey, A^ diffusing predators, and a circular haven of radius 
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R centered at the origin. Because of the recurrence of diffusion in two dimensions, 
the prey will eventually reach the haven if there are no predators, but the mean time 
to reach the haven is infinite. What happens when predators exist? How does the 
survival probability depend on the number of predators and on the initial positions of 
the prey and predators? How long does it take for the capture process to end? Another 
interesting two-dimensional geometry is a semi-infinite planar haven. Finally, in three 
dimensions, the transience of diffusion could lead to very different properties for the 
survival probability than in two dimensions. 
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